Hamiltonian approach for the wave packet dynamics: Beyond Gaussian wave functions 
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It is well known that the Gaussian wave packet dynamics can be written in terms of Hamilton 
equations in the extended phase space that is twice as large as in the corresponding classical system. 
We construct several generalizations of this approach that include non-Gausssian wave packets. 
These generalizations lead to the further extension of the phase space while retaining the Hamilton 
structure of the equations of motion. We compare the Gaussian dynamics with these non-Gaussian 
extensions for a particle with the quartic potential. 



Among various semi-classical approaches the Gaus- 
sian wave packet dynamics has been particularly suc- 
cessful and extensively used by physicists and chemists 
P, [3, H, 0, H, An especially attractive Hamiltonian 
formulation of this method was developed by a number 
of authors [1, S B S EH • In this formulation the Gaus- 
sian wave packet dynamics is recast into the classical-like 
dynamics that obeys Hamilton equations in the extended 
phase space that is twice as large as in the correspond- 
ing classical problem. The Hamiltonian approach offers 
a number of advantages over the non-Hamiltonian for- 
mulations [8| . Possibly the most important among them 
is that the existence of the classical-like Hamiltonian in- 
sures the stability of trajectories in the generalized phase 
space and well-defined quantum dynamics. A number of 
extensions of the Gaussian wave packet approach that 
include non-Gaussian wave functions either explicitly or 
implicitly were developed 11, 12, [T^j. Generally, how- 



ever, these extensions do not preserve the Hamiltonian 
structure of the equations of motion and, therefore, do 
not possess its advantages. The purpose of this paper 
is to develop a natural extension of the Gaussian wave 
packet dynamics that preserves the Hamiltonian form for 
the equations of motion. 

The Hamiltonian formulation for both Gaussian wave 
packet dynamics as well as its non-Gaussian extensions 
can the most easily be obtained through the time depen- 
dent variational principle. In this approach one intro- 
duces the functional 
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The requirement that ST = against independent varia- 
tions of (ip(t)\ and \ip(t)) leads to the Schrodinger equa- 
tion and its complex conjugate. Various approximation 
schemes can then be constructed by restricting \ip) 's to 
certain classes of functions that depend on a finite num- 
ber of time dependent parameters. In order to obtain 
Hamilton equations the parameters should be chosen in 
such a way that they satisfy the so-called canonicity con- 
ditions [141 • An alternative but related approach that 
does not introduces the canonicity conditions and leads 



to generalized Hamilton equations is developed in Ref. 

The derivation of the Hamilton equations for the 
Gaussian wave packets become very transparent if, rather 
then writing them in a specific representation, one uses 
their generic form as squeezed coherent states 0, [Toj | . 
Let us show how Hamilton equations are obtained for 
the Gaussian wave packets. These results will help us 
to demonstrate easily how to obtain the non-Gaussian 
generalizations . 

The Gaussian state is written in terms of the displace- 
ment and squeezing operators D(a) and S({3) as 



\G(a,P)) = D(a)S(P)\0), 



(2) 



where |0) is the ground state of the harmonic oscillator 
and displacement and squeezing operators are written in 
terms of creation and annihilation operators a and at as 



D(a) = exp(aa+ — a*a), S(/3) = cxp 
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(3) 

The state \G(a, (5)) depends on four real time-dependent 
parameters. These trial states are inserted in the func- 
tional ([T]) which will now be expressed in terms of these 
four parameters and their time derivatives. We then need 
to choose new parameters that satisfy the canonicity con- 
ditions. The choice of these new canonical parameters 
is completely determined by the differential part of the 
functional ([1]). Writing parameter a as r*iexp(i(/?i) and 
parameter (3 as r 2 exp^t^) and using the properties of a 
and aJ one can verify the following relations 
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Ds(a"<(e lipi coshr 2 - e ~^i-e2) sinnr2 ) 
-a (e"^ 1 coshr 2 - e 1 ^ 1 ^ 2 * 1 sinhr 2 )) , (4) 
iDs(rl + o+ri (e^ 1 coshr 2 + e -^fi-V2) silmr2 ) 
+ar 1 (e-' lipi coshr 2 + e 1 ^ 1 '^ sinhr 2 )) , (5) 
^S(a t2 e^-a 2 e-^), (6) 
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In the case of the Gaussian trial state ([2]) when calculat- 
ing the time derivatives for the approximating functional 
we need to multiply these expressions from the left by 
(DS) -1 and then average them over the vacuum state. 
Performing this averaging we obtain 



(G|A| G)=0 



(8) 



squeezed and displaced superposition of states with oc- 
cupation numbers that are multiples of 3. Recall that the 
(j^ squeezing operator acting on the vacuum produces states 
that are superpositions of number states with occupation 
numbers that are multiples of 2. In this sense the state 
\F) appears as a natural generalization of \G). We now 
have to choose some parametrization of coefficients c m 
that satisfies the canonicity conditions. Here we consider 
two possible types of parametrization for the function 

A simple two parameter parametrization of the state 
([T2")) that is close in spirit to the original Gaussian state 
is obtained by writing the state |£) as 



d 



(G\ — |G) = -sinhV (9) 

0(f2 2 



Using these results and introducing the new variables 
J\ = r\ and J2 = \ sinh 2 r2 while keeping the original 
tfiS we obtain the following form of the functional (fTJ) 
for the Gaussian state ([2]) 



J dti-^Jm-HdJihitpi})), 



(10) 



where H({ Ji}, {fi}) = (G\H\G) and h is set to unity. Pa- 
rameters Ji and ifi satisfy the canonicity conditions. The 
requirement that 5T = leads to the following Hamilton 
equations 
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dcpi ' 
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We now want to generalize these results to functions 
that are not necessarily Gaussian. An obvious general- 
ization of the squeezed state of Eq. ([2]) is the state of the 
form \F) = D(a)S((3)\£) where |£) is a state depending 
on an even number of real parameters. Note that when 
we consider the variation of functional T with function 
\F) instead of \G) we can still use Eqs. (jHEJ) for the 
derivatives with respect to r\ , ip\, r 2 and cf2 ■ To ensure 
that we end up with the Hamilton equations, we also 
want to preserve the property that the derivatives with 
respect to r% and r 2 that appear in Eqs. (|4|6|l vanish after 
averaging as in Eq. ((8|) . This implies that the state |£) 
must be such that the averages of a, a', a 2 , and <r are 
equal to zero. This is valid for any state that is a super- 
position of the number states with occupation numbers 
that differ at least by 3. In this paper we consider trial 
states |£) of the following general form 



M 
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that are assumed to be normalized with M being an in- 
teger that can range from 1 to infinity. Thus, instead 
of the squeezed and displaced vacuum as in the case 
of the Gaussian state, we are now dealing with some 
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where coefficients b m are fixed and the complex parame- 
ter 7 is allowed to change. The normalization constant N 
is given by JV( 7 ) = £~ |& m | 2 | 7 | 2m Eq. 03) is the gen- 
eral form for several classes of generalized squeezed and 
coherent states found in the literature, such as various 
multiphoton squeezed states [15( or eigenstates of oper- 



ator a d IpR Il8j . All these states are distinguished 
only by a particular choice of coefficients b m in Eq. (|13|) . 



Since these states are always normalized and reduce to 
the vacuum state for 7 = they can be written in terms 
of some two-parameter unitary operator T( 7 ) acting on 
the vacuum. The total trial state then has the form 
\F(a,/3)) = D(a)S(P)T(-y)\0) that looks like a direct gen- 
eralization of the squeezed coherent state D(a)S((3)\0) . 
The explicit form of T( 7 ) is only known for some spe- 
cial choices of 6 m 's [l5| and for this reason we will use 
the form (fT"3]) instead of T( 7 )|0) in the further analysis. 
Writing 7 as exp iip^ one easily verifies that 



{F\±\F) = 0, 

i 
{Fl d^J F) = Nil) 



(14) 
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= l -{F\(a^a)\F)=iK(r 3 ). (15) 

Using these relations as well as Eqs.(j4][7|) we obtain the 
new variables that satisfy the canonicity conditions and 
Hamilton equations 
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i(2X(r 3 )+l)sinh 2 r 2 , 



■h 



K(r 3 ), 
(16) 

with the original Pi's. The "classical" Hamiltonian 
H({Ji}, {<fi}) is given by (F\H\F). Application of these 
results to a specific Hamiltonian will require a particular 
choice for the form of coefficients b m in Eq. (fl~3"|) . Given a 
Hamiltonian, this choice can be guided by both physical 
plausibility and mathematical manageability. In particu- 
lar, we should be able to rewrite H. in terms of Ji 's rather 
then the original r^'s by using Eqs. (jTSJ) . Generally, this 
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implies that we should be able to solve the last of Eqs. 
([T6|) for r 3 . In the example below we take c m 's to be: 
c m = 1/Vmi. The relation of the corresponding state 
to the algebraic properties of the generalized three pho- 
ton creation and annihilation operators is discussed in 
Ref. [15J. We make this particular choice here primar- 
ily because of its mathematical convenience, since in this 
case K(r 3 ) = r 3 2 and, therefore, r 3 is easily expressed in 
terms of J3 . To distinguish this state from the other trial 
states we denote it 

Another, more detailed parametrization is achieved by 
treating coefficients c m in Eq. (fT2]) themselves as param- 
eters. To this end we write each of c m 's (for m > 0) 
as R m exp(i</> m ). By the normalization requirement cq is 

then given by cq = 
that 



1 - Ef l c " 



| 2 . It is easy to verify 
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(F\ — \F)=iRj 



(17) 



Using these relations as well as Eqs.(j4][7]) we obtain the 
new variables Ji that are expressed through the original 
parameters as follows 

00 

■h=rl J 2 = -(1 + 2^ 3m|i? m | 2 )sinh 2 r 2 , (18) 
and for i > 3 



m— 1 



Ji — Ri 



-2) 



•■Pi = 4>i 



(19) 



The original cpi and ip2 are the conjugate variables of 
Ji and J 2- The value of integer M determines by how 
much the phase space is extended compared to the Gaus- 
sian case. In the example below we consider the simple 
cases of M — 1 (six dimensional phase space) and M = 2 
(eight dimensional phase space). We will denote the cor- 
responding trial states by I-F2) and I-F3). An interesting 
open question that we do not consider here is how well 
can an arbitrary function be approximated by this type 
of the trial state when M is allowed to be infinite. 

Let us now compare these non-Gaussian extensions 
with the Gaussian wave packet dynamics in the case of 
the model Hamiltonian 
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(20) 



We will take the value of A = 1 and, therefore, the 
Hamiltonian is bounded from below. We will consider 
two choices for a: a = 1 and a = — 1. The latter case 
corresponds to a two- well potential with the two minima 
located at ±1. In order to make the comparison with the 
Gaussian approach we have to choose the initial state 
as a Gaussian one given by Eq. ((2]). We consider the 
same initial state for both choices of a with parameters 
a = and (3 — 0.1. This correspond to a slightly 

squeezed state with the average coordinate equal to one 
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FIG. 1: W(t) for \G) (gray dots and dashes), \F\) (solid gray), 
I.F2) (black dashes), and \Fz) (solid black) when a = — 1 . 



and zero average momentum. With these parameters the 
initial expectation values of the quadratic and quartic po- 
tential terms in Eq. (|20|) are of the same order. 

First we consider the case of a = — 1. The main quan- 
tity of interest for us is W(t), the squared overlap be- 
tween the exact wave function and its approximations 
given by \G),\Fi),\F 2 ), and \F 3 ) (Fig. The exact 

wave function is calculated by the numerical solution of 
the Schrodinger equation. To give the reader a better 
feeling of the time scales, Fig. [2] shows the time evolu- 
tion of the average coordinate (x(t)} over the same time 
interval. We can see from Fig. [T]that for short times (up 
to about 1 time unit) all three non-Gaussian extensions 
give better overlap with the exact wave function then the 
Gaussian approximation. For longer times, however, the 
behavior of W{t) for different approximations becomes 
more complicated. Depending on a specific time each 
of the approximating wave functions (with the exception 
of I-F2)) can give the best overlap with the exact wave 
function. When averaged over the time interval of the 
Figure (20 time units) the values of the squared overlap 
are the following: 0.466 for \G), 0.612 for |Fi), 0.450 for 
IF2), and 0.492 for \F 3 ). Thus, on average for this time 
interval |.Fi) provides a much better description then the 
Gaussian wave function, \F 3 ) performs slightly better, 
and |_F2) gives worse results then the Gaussian. These 
results are in agreement with Fig. ([2]) which shows that 
l-Fi) generally reproduces exact (x(t)) better then the 
other functions. 

The squared overlaps for the case of a = 1 is shown 
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FIG. 2: (x(t)) for the two-well potential (a = —1) for the 
exact wave function (long gray dashes), \G) (gray dots and 
dashes), \Fi) (solid gray), IF2) (short black dashes), and IF3) 
(solid black). 

on Figure [3] We see again that for up to about 1 time 
unit the non-Gaussian wave functions are closer to the 
exact one then the Gaussian wave packet. For longer 
times no trial function seems to offer a definite advan- 
tage. W(t) averaged over 20 time units has the following 
values: 0.353 for |G), 0.353 for |iq), 0.336 for \F 2 ), and 
0.338 for \F 3 ). For this case both \F 2 ) and \F 3 ) give 
slightly worse averaged W{t) then the Gaussian case. 

To sum up, in both cases (a = ±1) the short time 
dynamics is improved by using extended wave functions. 
For longer times \F\) performs at least as good as \G) 
or better, while \Fi) and \F%) show no definite improve- 
ment over \G) and can even perform worse then \G). This 
long time behavior seems surprising. Indeed, one usually 
expects that variational approaches give results that be- 
come progressively better when the number of variational 
parameters increases. We do not offer an explanation of 
this kind of behavior in this paper. It is possible that the 
reason for this should be sought in the interpretation of 
the variational principle when applied to the functionals 
that linearly depend on the time derivative of the field. 
When performing the variation of the functional |T|) to 
derive the Schrodinger equation we assume that we can 
independently vary \tp) at both time limits. In the re- 
sulting Schrodinger equation, however, the initial wave 
function uniquely determines the final one. 

To conclude, we constructed several generalizations of 
the Gaussian wave packet dynamics that include non- 
Gaussian wave functions and can be formulated in terms 



FIG. 3: W(t) for \G) (gray dots and dashes), \Fi) (solid gray), 
IFb) (black dashes), and [JF3) (solid black) when a = 1. 



of the classical dynamics in the extended phase space. 
Applied to a model system these approaches give a bet- 
ter description of the exact wave function for short times. 
Understanding the peculiarities of the longer time behav- 
ior as well as further extensions of these approximations 
require additional research. 
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